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Abstract 

We compute the spatial-volume dependence of the spectrum of 4D SU(3 < < 6) gauge 
theories by lattice Monte-Carlo techniques. Setting the scale with the string tension, the 
spatial volume is Li^ with 0.78fm < L < 2.3fm. The Euclidean 'time' direction is kept 
large enough to be considered infinite and the boundary conditions are periodic in all four 
dimensions. We study the mixing of torelon pairs with the scalar and tensor glueballs, using a 
2x2 Hamiltonian based on large- counting rules. Looking to the other symmetry channels, 
finite- volume effects on the glueball spectrum are already surprisingly small in SU(3), and they 
become rapidly smaller as is increased: several low-lying SU(6) states have no finite-volume 
corrections at the 1-2% level, at least down to L = 0.9fm. We discuss the relation of this 
work with analytic calculations in small and intermediate volume, and with Eguchi-Kawai 
reduction in the planar limit. 



^On leave from: Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford 0X1 3NP, U.K. 



1 Introduction 



It has long been realised that the volume dependence of SU(A^) gauge theories is an important 
source of information on their dynamics. In very small spatial volumes LA <C 1, perturbation 
theory is applicable due to asymptotic freedom and dimensional reduction effectively occurs 
with respect to all three space dimensions (we assume periodic boundary conditions). This idea 
was used in to derive an effective Hamiltonian for the zero modes of the pure gauge theory, 
which allows one to compute the 'glueball' spectrum analytically. The technique was then 
gradually extended by taking into account tunneling between perturbative vacua to reach 
volumes corresponding roughly to 0.5fm [2| (for a review, see After pioneering small- 

volume Monte-Carlo calculations by Berg and Billoire the analytic results were found 
to be in agreement with numerical data for the SU(2) gauge group in the much-improved 
calculations El of the late eighties; numerical results for the SU(3) gauge group remained 
more limited [7j. Also in the presence of light quarks, finite- volume techniques have proven 
useful to gain insight into the theory. An effective small- volume action was derived in |^ for 
the SU(2) gauge theory including the one-loop contribution of light quarks, and later extended 
to intermediate volumes 0. The existence of analytic results in small volume, even at finite 
lattice spacing was exploited in ^01 to study the scaling behaviour of various lattice gauge 
actions. More recently, the so-called e-regime, defined by ni.,^L <C 1, has been used as a testing 
ground for chiral perturbation theory, through the connexion with random matrix theory . 

At intermediate volumes, L ~ Ifm, non-perturbative effects are an obstacle for purely 
analytic calculations. However, physical quantities of the finite-volume theory can be related 
to observables of the corresponding infinite-volume theory. The finite-volume correction to 

the mass of a stable state is proportional to e~~"^^ (where m is the mass gap) and the 
proportionality constant is nothing but the forward scattering amplitude of two of the lightest 
particles in infinite volume jl2 | ll,Sj. Further, unstable states will mix with the scattering states 
of their decay products in a finite volume, and the mass splittings inform about the phase 
shifts in infinite volume, and in particular on the width of the resonance |14( I15j. Similarly 
matrix elements in finite volume can be related to the 'weak' decay rate of 'kaons' in infinite 
volume '16'. 

Another point of view on gauge theories in finite-volume is obtained from t'Hooft's large- 
expansion, where is the number of colours Indeed, large- A^ QCD is thought to 

be a theory of free mesons and glueballs (while baryons have masses of order N |18j): this 
picture is supported by arguments involving the large- A^ counting rules for Feynman diagrams 
and the assumption of confinement at A^ = oo. From the general relations between finite 
and infinite observables cited above, one immediately infers that finite-volume effects should 
become weaker as A^ increases, since widths and scattering amplitudes are suppressed by a 
power of 1/A^. Conversely, while the low-lying SU(A^) glueball spectrum in large volume has 
been shown to have 1/A^^ corrections . the finite- volume effects on energy levels can 

be used to determine non-perturbatively the A^-dependence of certain three- and four-point 
functions. This work can be seen as a first attempt in this direction. By comparison, measuring 
these n-point functions on the lattice in a straightforward way would be extremely hard. 

The argument on the size of finite-volume effects outlined above fits in with the ideas of 
reduction initiated by Eguchi and Kawai the SU(oo) lattice gauge theory defined on a 
single lattice hypercube of size has equivalent loop equations for its Wilson lines to those 
of the corresponding theory on an infinite lattice, and corrections come from the connected 
correlation of traced Wilson loops, which is suppressed by 1/N'^. In fact the proof given 
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by the authors trivially extends to any lattice, L = La. Thus if the vacua of the theory, 
defined (a) on the lattice and (b) on the infinite lattice, 'choose' the same class of solutions, 
finite- volume effects on Wilson loop expectation values vanish at N = oo. 

The equivalence of the loop equations holds provided that the centre symmetry Z{N) 
(which is promoted to a global U{\) symmetry at = oo) does not break spontaneously. 
Strictly speaking, at finite such a spontaneous symmetry breaking cannot occur in finite 
volume. But the energy barrier between the different vacua allowed by the symmetry grows 
with N so that tunnelling between them becomes completely suppressed at = oo. This point 
can be made more precise. In the Euclidean path integral formulation one may reinterpret one 
of the 'small' spatial dimensions as the Euclidean 'time' direction, which fixes the temperature, 
T = 1/L > Tc- In this picture, two of the three spatial dimensions have size L and one is 
infinite (or at least > L). The thermal tunnelling from one vacuum to the other was discussed 
in in a semi-classical way. In this process two different vacua must coexist 'temporarily', 
and a domain wall must appear between the confined and one of the N deconfined vacua. 
Once the 'temperature' 1/L is large, perturbation theory becomes applicable. The domain 
wall tension between two Z{N) vacua is given by |21j 

If a domain wall orthogonal to the large dimension separates two A;- vacua, its energy is a^L"^ ■ 
At temperature T = j^, the tunnelling rate is thus proportional to the Boltzmann factor 
g-bk(N,L) ^ We note that due to asymptotic freedom, even at fixed A^ the Boltzmann factor 
is suppressed as L ^ 0, albeit more slowly than a power law. To obtain the tunnelling rate 
this expression ought to be multiplied by further factors that take into account the entropy 
of domain walls; however they are not expected to depend strongly on A^, and therefore we 
conclude that the tunnelling is exponentially suppressed in A^. At large A^ the exponent is 
large, and this justifies the semi-classical method and language used in the argument. If we now 
return to the original interpretation of the system as a finite- volume one at zero temperature, 
we see that from the point of view of Monte-Carlo simulations, the centre symmetry along 
one of the dimensions of size L is practically broken permanently at finite, large A^. This 
conclusion does not depend on the size of the dimension left unspecified, in particular it holds 
for and L'^ x oo lattices. 

It is well-known that if one takes L = 1 the spontaneous symmetry breaking does oc- 
cur (in the sense outlined above) with periodic boundary conditions, a fact that led to the 
'quenched' j22j and 'twisted' Eguchi-Kawai models. In the former case, the path integral 
is constrained in a way that maintains a uniform distribution of eigenvalues of the link ma- 
trices on the unit circle. In the latter case, the symmetry breaking is prevented by twisted 
boundary conditions |2ZI- In fact, the connexion between the twisted reduced model and the 
ordinary gauge field theory is then even more precise: for the choice of symmetric twist, the 
reduced theory at finite- A'' has loop equations equivalent (up to the 0(1/A''^) corrections to 
factorisation) to those of the corresponding theory on a periodic lattice, with = y/N. 
In particular, at large but fixed N, (TrC/^) = and {{TrU^)^'^) = 0, but the latter (a Polyakov 
loop in the field theory) will have a very long autocorrelation time at sufficiently weak cou- 
pling. Therefore, as one takes the continuum limit, A^ has to be increased so as to maintain 
the physical size of the box L = aL ^ in which the gauge theory is defined in order to converge 
to the expected theory. 
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To compare the spectrum of the (partially) reduced and the original theory, connected two- 
point functions must be considered. A Hamiltonian version of the reduction (on an x R 
lattice) was discussed in |26j and arguments for the equality of the spectra in the reduced 
model and in the infinite-volume theory were given. In essence, one of the arguments is 
that the reduction property still holds if a source term for closed Wilson loops is added to 
the action; since this is a generating functional for glueball correlation functions, the finite- 
volume corrections to the spectrum must also vanish in the planar limit. For the argument 
to be complete, however, the equality of the generating functionals must be shown to hold 
beyond the linear term in the source. 

Factorisation implies that the N = oo theory corresponds to a classical theory. In the 
coherent state method [211, the spectrum could in principle be obtained by expanding the 
classical Hamiltonian (including the 0(l/iV^) corrections) to quadratic order around the point 
of phase space corresponding to the vacuum state. 

In this work we investigate by Monte-Carlo techniques the finite-spatial-volume spectrum 
of the SU(3), SU(4) and SU(6) gauge theories, for 0.78fm <L< 2fm. By keeping L > 1/Tc, 
we will avoid that the system gets locked into a Z{N) broken phase for long periods of Monte- 
Carlo time (especially for 5[/(4) and SU{Q)). We study separately the symmetry sectors 
where 'additional' states occur and those where there are only 'intrinsic' finite-volume effects 
on the states already present in the infinite-volume theory. In particular we examine the 
dependence of the finite-volume effects in both cases. 

In the next section we describe in more detail the theoretical expectations on the behaviour 
of the finite- volume spectrum of pure gauge theories. We present the lattice results for the 
SU(3) gauge group (section OJ, and then SU(4) and SU(6) (section 0}. We finish with a 
summary of the results and a general discussion. 

2 Finite-volume and planar physics 

We consider a four-dimensional Euclidean quantum field theory defined on a spatial hyper- 
torus with an infinite time direction, y = L'^ x R. In a quantum field theory admitting a 
mass gap m, the leading finite-volume effect on the spectrum normally corresponds to the 
emission, propagation 'around the world' and reabsorption of the lightest particle. A formula 
corresponding to this picture was established to all orders in perturbation theory in |12|: the 

mass shifts behave as e~ 2 and the formula is accurate when the exponent is large. A 
simple picture for the origin of this mass shift is that the particle interacts via a Yukawa 
potential with its 'replica' located by periodicity in L x L x L cells. It is now tempting to 
jump to the conclusion that the spectrum only has exponentially small corrections. However 
the situation changes qualitatively if a 'new' light state appears at volumes where L is still 
large compared to the inverse mass gap of the infinite-volume theory - a non-perturbative 
effect. This is precisely the situation of pure SU(A^) gauge theories, where such an extra state 
is present in a finite box but becomes a very massive and broad resonance which decouples in 
the infinite-volume limit. We now review in some detail how this happens. 

Due to the global Z{N) symmetry in the SU(A) gauge theory, the eigenstates of the 
transfer matrix can be classified according to their winding number. Let us encode the Z{N) 
transformation properties of an operator or state by a triplet {ex-,ey,ez) ^Sl- A glueball has 
(0,0,0), a torelon in the x direction has (1,0,0), the triplets combine additively for a direct 
product of states and the entries are defined modulo A^; moreover, charge conjugation connects 
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the sector k to N — k. In particular, it turns out that the hghtest state for L < 1.5fm is the 
torelon, which acquires the phase e*27r/7V ^j^jgj- ^ centre transformation. One can then imagine 
a state made of two torelons winding the same cycle of the hypertorus in opposite directions. 
Such an object is invariant under the Z{N) symmetry, and therefore becomes part of the 
glueball spectrum. The three resulting states (one for each spatial direction) can then be 
classified into the singlet Af~^ and the doublet E^^ of the cubic irreducible representations. 
In general there is a splitting between these two 'extra glueball states'. 

Apart from the simple torelon pairs, one could imagine that even more states appear in 
the (0,0,0) sector in small volumes. For SU(iV), N torelons winding a cycle of the torus in 
the same direction transform trivially under the centre of the gauge group. The real part of 
such operators couple to the and sectors, while the imaginary part belongs to the 

triplet. For k < N/2, two {k, 0, 0) operators constitute another possibility to construct 
glueball operators. One can also consider operators transforming non-trivially under two 
or three of the Z{N) symmetries, (1,1,0) and (1,1,1); two of them, winding in opposite 
direction, again couple to (0, 0, 0) states. However, in the absence of interactions these states 
are all significantly heavier than the lightest two states in the representations we shall consider. 
In the following we shall see no numerical evidence that they play a role in the dynamics of 
the low-lying spectrum. 

At L ~ 0.8fm, the single-torelon mass is significantly lighter than half the lightest scalar 
glueball mass in large volume. Therefore, in the absence of interactions, the lightest state 
in the Af^ (and E^^) channels would be a two-torelon scattering state. The interaction 
potential between the two zero-momentum states is in general suppressed by their dilution in 
the transverse directions, AV oc {£±/L)'^ in our case, where £± is the average width of the 
'string'; but this is hardly a suppression at the volumes we are interested in. Also, (working 
in the basis of states with definite relative momentum) the higher relative momenta will 
contribute little to the ground state for L < Ifm. 

Since the torelon mass varies roughly proportionally to L and the glueball mass is only 
weakly dependent on L, we expect the dynamics of torelons and their mixing with an ordinary 
glueball state to determine the mass of the two lightest Af~^ (and £'++) states in a certain 
range of volumes. It turns out that it corresponds roughly to 0.8fm < L < 1.2fm. In that 
range, there is a large energy gap to the third state in these two representations. That means 
that to a good (and controllable) approximation, we can study the system with a 2 x 2 effective 
Hamiltonian. 

In order to apply the stationary perturbation theory of quantum mechanics, we need to 
choose a starting orthonormal basis which diagonalises an 'unperturbed' Hamiltonian. We 
suppose here that the full Hamiltonian of the SU{N) gauge theory can be expanded in inverse 
powers of l/AT: 

HiL,N) = ±^. (2) 

k=0 

The existence of the t'Hooft limit implies that Hq has the same eigenvalues as the Hamiltonian 
of the SU(oo) theory^; the latter, according to conventional wisdom, describes non-interacting 
colour singlets. Ho{L) is our 'unperturbed' Hamiltonian. If mT{L, N) is the torelon mass and 
mG{L,N) the glueball mass in the relevant symmetry channel, then the restriction of Ho{L) 

^The SU{oo) theory may have additional, decoupled sectors of free hadrons. 
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to the vector space span by these two states reads 



H.{L) = (^'^^f°'^ ° ). (3) 

\ mciL, oo) J 

We now consider the fuh hamiltonian H{L, N) in the same basis: 

rr (T ^^_( E2t{L,N) h{L,N) \ 
H^x^(^'^)-[ hiL,N) maiL,N) ) 

with 

E2t{L, N) = 2mT{L, oo) +vt{L,N), rndL, N) = mciL, oo) + vg{L,N). (5) 

The diagonal matrix elements vt and vq contain the terms describing 'intrinsic' finite-volume 
corrections of the unperturbed states. For instance, in an effective low-energy description, the 
force between flux-tubes is attractive as a result of the exchange of the lightest particle, a 
scalar glueball; such a force is always attractive. The glueball state also undergoes 'intrinsic' 
corrections, proportional in leading approximation to e~'^"^^ However, equally important 
are the off-diagonal elements describing mixing. Pictorially one can imagine the following 
process: two flux-tubes undergo 'fusion', after which they form a closed, contractible flux 
loop. In the flux-tube model j^H], a glueball precisely corresponds to such a loop: at this point 
the uncertainty principle (and possibly curvature effects) prevent the loop from shrinking to 
a point. 

Let us now discuss the N dependence of the various contributions to the Hamiltonian Q). 
The finite- energy-shift vq corresponds to the usual correction to the mass of a colour- 
singlet state, along with the virtual emission 'around the world' and reabsorption of a colour- 
singlet by a colour-singlet state, which is also a effect. Similarly vt{L) contains 0{l/N'^) 
corrections to the mass of each torelon, and the 'elastic' interaction of torelons encoded in 
vt{L) proceeds through the exchange of a colour-singlet, and therefore is also 0{1/N'^). h{L) 
on the other hand, corresponds to the decay amplitude of a colour-singlet into two colour- 
singlets, which is 0{1/N). A simple way to see this is that the decay rate of a glueball above 
the two-torelon threshold (if the transverse dimensions were large) would be proportional to 
h?' in this formalism, according to 'Fermi's golden rule'. On the other hand, the width of 
colour-singlets in the pure gauge theory is 0(l/iV^), which leads to the announced conclusion. 

Thus we expect, at sufficiently large A^, the functions vg, vt and h to take the functional 
forms 

vg{L,N) ~ -^vg{L) (6) 

vt{L,N) ~ -^vt{L) (7) 

h{L,N) ~ ^HL). (8) 

On the lattice, we of course compute what corresponds to the two eigenvalues of H2x2iL, N), 
mo{L,N) and mi{L,N). Suppressing everywhere the {L,N) dependence, one has 

mo + mi = E2T + rriG (9) 
{nio-rmf = {E2T - mcf + 'ih^ (10) 
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The information extracted from the lattice is in general not sufficient to determine all three 
unknown functions vq, vt and h. However the large- counting rules can help us out. 
Firstly, according to the latter, 

mT{L,N) = mT{L,oo)+0{l/N^). (11) 

For the glueball mass, similarly 

maiL = oo,N)= mciL = oo, iV = oo) + 0{l/N^), (12) 

since in large volume there is no other light state the glueball can mix with, and 

mG{L = oo, N = oo) = mciL, N = oo) (13) 

if the planar theory has no finite- volume corrections. Therefore, using Eos. 1511^ 171 II 11 IT^ and 
[T3l Eq. ini becomes the prediction 

mo(L, N) + mi(L, N) = 2mT{L, N) + mciL = oo,N) + 0{l/N^) (14) 

that can easily be tested on the lattice, since all quantities are at a given N. 

Secondly, the half-splitting of the two lattice states provides an upper bound for h{L, N): 

\h{L,N)\ <^\mo{L,N)-7m{L,N)\. (15) 

For L such that 

Eq. ^1 tells us that h{L, N) saturates the upper-bound, up to small corrections. Since 
E2t{L, N) and mciL, N) are not known a priori, it is more convenient to express the condition 
above by trading them for the known quantities 2mT{L, N) and mciL = oo,N): 

1 - ^^^#4t « Od/iV), (IT) 

mcioo, N) 

(the left-hand side of [H and HZI only differ by 0{l/N^) terms). If h{L,N) is found to be 
0{1/N), the procedure used to extract h is self-consistent. What is interesting, from the 
point of view of 1/A^ counting, is that in the volumes satisfying Eqn. El the finite- volume 
mass shift is of order 1/A^. Because in large volume glueballs are relatively heavy in units of 
the string tension, it is easy to see that this condition will necessarily be satisfied at some 
intermediate box size. Thus there are circumstances where the expectation that finite-volume 
mass shifts are suppressed by can be wrong. Analogous violations of the '1/iV^ correction 
rule' for the spectrum of SU(A^) gauge theory will be discussed in section |SJ 

In conclusion our expectations are the following. In most representations of the cubic group, 
where no degeneracies appear as the volume is varied, the energy levels have small deviations 
from their infinite- volume values. These corrections are of order and behave as e~'^™(^)^ 

(as long as the exponent is large). In those representations where degeneracies appear, such as 
the Af^ and which contain torelon pairs, an interesting interplay between the parameters 
L and A'^ appears. When L is such that the relative separation of states is much larger than 
0{1/N), the deviations of the glueball masses from their infinite- volume values are 0(l/iV^). 
But as the torelon mass comes within 0{1/N) of half the infinite- volume glueball mass, the 
finite- volume corrections to glueball masses are enhanced through mixing to 0{1/N). In 
parallel, as the torelon-pair energy falls below the lightest ^4^"*" glueball, the mass gap m{L) 
is given by that energy and finite-volume shifts become more rapid, as a function of L. 
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3 Lattice results for SU(3) 



Our lattice calculations employ the standard plaquette action and the equally standard 1:4 
combination of heatbath and over-relaxation sweeps is used for the update. We calculate 
ground and excited state masses m in the zero-flux sector (0, 0, 0) from Euclidean correlation 
functions of zero-momentum operators using standard variational techniques jSlj . We extract 
the string tension, a, from the mass of a flux loop that closes around a spatial torus. We 
perform calculations for a fixed value of the inverse bare coupling /? = Q/g'^ = 6.0 corresponding 
to a lattice spacing a ~ 0.097fm, if we use ^/a = 440MeV. The calculations are on lattices 
ranging from 8^48 to 24'^24, corresponding to a spatial extent 0.78fm < L < 2.3fm. Since some 
states become very light on small volumes, the time direction is gradually increased from 24 
to 48 lattice spacings as the spatial volume is reduced in order to suppress 'around the world' 
contributions (proportional to e~™'(^*~*)) to the correlation function C{t). In the L = 8 case 
we checked that increasing Lt to 64 does not alter the measured spectrum. 

It is well-known that it is crucial to include sufficiently different operators in the variational 
basis that overlap significantly onto the low-lying spectrum, particularly when one is trying 
to extract more than one energy level. From the discussion of the previous section, it should 
be clear that what is required in the Af^ and E~^^ sectors are contractible Wilson loops as 
well as pairs of spatial Polyakov loops winding around a cycle of the spatial torus in opposite 
directions. We measure them with zero relative momentum, since a crude estimate of the 
energy of two non-interacting torelons with one unit of momentum 2tt/L each turns out to 
be well above the lightest Af~^ state for the L of interest^. We take the linear combinations 
of the scattering states that belong to the Af^ and E~^~^ representations of Oh- We note 
that the overlaps between the two-torelon and the Wilson loop operators is substantial in 
the region L ~ Ifm (see below). In the other representations, our basis of operators consists 
only of appropriate linear combinations of contractible Wilson loops. All our operators, the 
Wilson loops and the scattering state of Polyakov loops, are constructed with several levels of 
smearing and blocking (see e.g. ^30J. We note however that one must be careful in a small 
volume not to include operators so non-local that they contain paths winding around a cycle 
of the torus. 

In some of the simulations we use a 2-level algorithm [H2j that reduces the variance of 
rapidly decaying correlators, thereby allowing us to extract the mass at larger time-separation. 
This is especially relevant in small volumes, where it turns out that our fuzzy operators do not 
necessarily have as good overlaps on the physical states as in large volumes. On the smallest 
volume we notice a significant increase in the auto-correlation time of the spatial Polyakov 
loop expectation value (see Fig. [2 both of these runs were employing the ordinary algorithm). 
This may be attributed to the proximity of L to 1/Tc. aTc ~ 7.25 according to [20], and it 
appears that tunnelling transitions are moderately frequent on the 8^ x 48 lattice. In such a 
situation a possible criterion to split up the data into independent blocks is to require that 
the magnitude of the average Polyakov loop in each block be less than a prescribed value. 
Our jackknife bins cover a Monte-Carlo time that is larger than the width of a typical peak 
on Fig. [21 The way we store the measurements J22| in the 2-level algorithm allows us to study 
in detail the dependence of the statistical errors on the jackknife bins; we are confident that 
our statistical errors are not biased by more than 10%. 

Table ^ contains our data on the various cubic representations. In the next section we 

■^We reserve the name 'adjoint Polyakov loop' for the case where the Polyakov loops are created at the same 
position. The use of such operators was pioneered by Berg and Billoire 
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study the lightest few states in the Af^ ('scalar') and E^~^ ('tensor') channels before we move 
on to the other representations. 



Channels with the lightest states 

The evolution of the lightest Ai~^ and E~^~^ states is illustrated on Fig. El Although we are 
only really interested in the first two states, we monitor the level of the third excited state to 
estimate the mass gap to the next states. We also plot the T^~^ states which are degenerate 
with the E~^^ states in the large volume (and continuum) limit. 

We observe that below 1.2fm, there are two states under 2GeV, rather than one; they 
are the Af^ and ^4^"*"* states. Meanwhile the Af^** is approximately degenerate with the 
large- voume Af~^*. 

The straightforward interpretation, namely that the fundamental state becomes lighter, 
and that the n*^ excited state becomes degenerate with the large- volume (n — 1)**^ excited state 
(n = 1, 2, . . .), would be hard to understand in a simple way. A far more natural interpretation 
is the one developed in the previous section, namely that one particular state, associated with 
the torelon pair, becomes light on the small volumes. It crosses the level of the large-volume 
Af~^* state around L = 1.5fm, and that of the large-volume Af~^ level around L = Ifm. Each 
level-crossing is accompanied by a maximal mixing of the torelon pair with the glueball state. 

A fit to the lightest energy level for L > 12 was attempted, assuming the functional form 
to m{L) = moo — X ^-^P (~^"^oo-^^) |l2j' The fit is acceptable but we find that the factor C 

is large (670 it 150) and driven by the point at L = 12, signalling that the glueball exchange 
'around-the-world' is probably not the dominant effect in that range of volumes (such a large 
factor was already found in [HSl)- The two lightest states come closest together at L = Ifm. In 
a two-state mixing model, the identity of the two states at smaller L is inverted: the ordinary 
glueball component of the first-excited state is larger than that of the ground state, which is 
now more of the torelon-pair nature. This crossing apparently happens when the curve mo{L) 
has an inflexion point. This interpretation is supported by the fact that our two-torelon 
operators have significant overlaps^ with the usual scalar glueball operators (0.32 at L = 11, 
0.67 at L = 10 and 0.55 at L = 9, choosing each time our best operator of each type). 

A similar pattern is seen in the channel: the ground state is pushed down by the 

presence of an extra light state for L < 1.2fm. Below that L, the first excited state is roughly 
degenerate with the fundamental T^"*" state, which has only small finite-volume corrections. 
The same mixing and crossing interpretation thus seems appropriate: on the smaller volumes 
it is the first-excited state which has a glueball-dominated wave function. 



Two-state mixing model 

At L = 12, we observe a state with an energy of almost exactly 2mT{L), both in the Af^ 
and E~^^ sectors, and the other light state still has very small finite-volume corrections with 
respect to the large-volume glueball masses. However the situation changes quite rapidly 
for L < 11. Compared to other symmetry channels, the scalar and tensor glueballs masses 
suddenly undergo finite-size shifts larger by about a factor of 2. 

The near-degeneracy of the two lightest states in the region 0.9 — l.lfm suggests the appli- 
cation of the two-state mixing model. Fig. ^ shows the average mass of the fundamental and 
first excited states (separately in the Af~^ and sectors). The lattice data is compared 



The overlap is defined as 



(O1IO2) 
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to the energy level mT{L) + 1710(00) /2 that one expects if torelon interactions and intrinsic 
finite- volume corrections to the glueball mass vanish (see Eq. I14() . The agreement is good 
for L > 0.9fm, while at smaller L the true average mass is pushed up compared to the naive 
estimate. This provides evidence that the 'intrinsic' effects are small in the region of interest, 
though one could be fooled if the mass shift for the torelon-pair has opposite sign to that of 
the finite-volume glueball mass, which would result in a large cancellation of these effects in 
the quantity considered. 

Ignoring this possibility for the moment, if we now suppose (in the spirit of large- count- 
ing rules) that the 'unperturbed' energy levels are separately well described by the values 
2mT{L) and mG(oo), we can estimate the mixing strength h[L) through Eqn. Ilfll For L = 10 
we find h{L) = /i ~ 0.15(2)mG'(oo) (the quoted error is statistical). In the E'^^ sector, the 
splitting is smaller than ImG'(oo) — 2mT{L)\ in the 0.9 — l.lfm region, and mixing can only 
split the states further apart. Perhaps unsurprisingly we must conclude that = 3 is not 
yet large enough for the A^-counting rules of section 2 to be reliable in the extraction of the 
mixing strength h. Nevertheless the mass splitting, given in Table^l provides an upper bound 
on h{L). 

Other representations 

In the other representations the complications discussed above do not occur. We consider 
the r2++, E-+, T^^, A--^ and representations (see Figs. El EJ. The T^^ state is the 
easiest case, because of its relatively small mass. As mentioned above, its volume dependence 
is weak and smooth. The same holds for the E '", T2~^ and states, whose masses increase 
monotonically as the volume becomes smaller. The increase is at the 10% level at L ~ 0.88fm, 
as we discuss more quantitatively below. Unfortunately our calculation in the 'pseudoscalar' 
channel A^^ is not sufficiently good to allow us to make precise statements about its volume 
dependence, but it is clear that it is the only one of these states that is becoming lighter - its 
mass seems to fall off rather abruptly around Ifm. We remark that our operators had poor 
overlap onto the state on the smaller volumes, whilst having satisfactory overlap on the large 
volumes. 

4 Larger N 

We repeat the calculation of the spectrum for SU(4) and SU(6); our method is entirely anal- 
ogous to the SU(3) case. The lattice spacing is slightly larger, a = 0.1 Ifm if we set the scale 
with the string tension, but prior knowledge |33j indicates that the system is already in the 
scaling region at that lattice spacing. In the Af~^ and E^~^ sectors, we again include the 
two-torelon operators along with the Wilson loops. Zero- momentum scattering states of two 
loops with AA-ality k = 2 were not included, although they could have a substantial overlap 
with one of the first excited scalar glueballs for A^ > 4, if its wavefunction were well described 
by a closed k = 2 string. This possibility has recently been suggested but we leave it as 
a subject of investigation for the future. The SU(4) and SU(6) spectra are given in Table 121 
and 131 

A^-dependence of 'intrinsic' finite-volume effects 

We start with the simpler case, namely the comparison of states not affected by the torelon 
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pairs. On Fig. [3 (top) we compare the relative finite- volume mass shifts 



of the lightest state between L ~ 0.85fm and L = 1.2fm for the representations Tg"^, , 
E ^ and {Loo corresponds to our largest lattice). Note that the error bars on m{L) and 
m{L = oo) are independent and therefore combine in quadrature. We already observed in the 
previous section that the variations were small in SU(3), at the 10% level for L ~ 0.88fm. 
Here we see that the finite-volume dependence gets even more suppressed as one moves from 
SU(3) to SU(4) and SU(6). This suppression is rapid: for = 6 our data is consistent with 
vanishing finite- volume effects, in spite of relatively small error bars (1—2% for the T^"^ and 
cases). Thus determining the N functional dependence numerically will require very 
accurate data. On the other hand the finite-volume corrections vary quite smoothly with L 
in these representations, and therefore it is sufficient to compare different SU(A^) theories at 
one chosen value of y^L ~ 2, for instance. 



dependence of torelon-pair mixing with glueballs 

As we showed in the SU(3) case, the spectrum in the A^^^ and E^^ sectors is affected by 
the presence of torelon pairs. The average mass and the mass splitting, that appear on the 
left-hand-side of Eqs. Inland are given in Table 4. Recall that the half-splitting provides an 
upper-bound for h{L,N) (Eqn. [T5|) . Our data shows that the mixing strength is modest in 
magnitude but the determination of its N dependence requires accurate data for several L in 
the region of maximal mixing given by Eqn. 1171 Note that the minimum of the half-splitting 
as a function of L gives h only if the latter is L independent. 

Considering the bottom part of Fig. [3 we see that the finite- volume mass shifts, although 
smaller than for SU(3), remain significant also for SU(4) and SU(6) in these representations, 
and that they are strongly dependent on the precise value of L. Due to the crossing of states, 
the shift changes sign around Ifm (we used the state which is closest to mciL = oo) to 
compute the relative mass shift). These observations are consistent with the expectations on 
the dependence of finite- volume effects formulated in section 2. 



5 Discussion 



We have studied the finite-size evolution of the low-lying glueball spectrum as a function of the 
number of colours in SU(A^) gauge theories. The box is in size, with 0.78fm < L < 2.3fm, 
the boundary conditions are periodic and the lattice spacing is fixed at roughly O.lfm. Seven 
symmetry channels were investigated. In most of them (Tj"^, E T2~^, ), the states 
observed connect smoothly with their infinite- volume counterparts. In the SU(3) case they 
become gradually more massive, and the mass shift is only at the 10% level at 0.8fm. For 
SU(4) and SU(6) this mass shift is very suppressed, and the SU(6) data is in fact consistent 
with vanishing finite- volume effects. This constitutes the first observation of the large- iV 
suppression of these effects on the spectrum in Monte-Carlo simulations. To show that these 
effects are truly of order 1/A^^ will require very accurate data. We also note that asymptotically 
Liischer's formula jl2| I13j predicts a lowering of the glueball masses as L decreases, since the 
exchange of the lightest scalar glueball yields an attractive force; we conclude that the effect 
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we are seeing is a different one, and that the asymptotic correction is too smah to be observed 
with the present statistical accuracy. 

In the A'^^ 'pseudoscalar' sector, the hghtest glueball becomes significantly lighter with 
respect to the infinite- volume state, and since this is observed for all gauge groups considered 
we think that this state deserves a dedicated study. One might speculate that the mass shift 
is related to the lower average in small volumes {Q is the topological charge). In that 
case, one would also expect the pseudoscalar glueball to show up as a lighter state in the 
presence of light dynamical quarks, which suppress the topological charge fluctuations. We 
note that a surprisingly sharp onset of the validity of random matrix theory in predicting 
the distribution of eigenvalues of the overlap Dirac operator in different topological sectors 
on quenched configurations was observed around l.lfm in It is not excluded that the 

behaviour of the lightest pseudoscalar glueball is related to that observation. 

The Af-^ (scalar) and (t ensor) sectors are most simply interpreted in terms of the 

appearance of extra light states, with respect to the large volume situation. These have a 
natural interpretation in terms of pairs of flux loops winding around one cycle of the hypertorus 
in opposite directions. Around Ifm the states come close to the ground states in these sectors. 
A simple two-state mixing model was used to interpret the two lowest energy levels. It appears 
that the binding energy of the torelons and the intrinsic mass shifts of the glueballs are rather 
small, so that the main cause of the variation of the spectrum in the L = Ifm region is mixing. 
In SU(3) it was estimated that the mixing energy amounts to about 15% of the infinite- volume 
glueball mass in the scalar sector. This estimate still has a very large systematic uncertainty; 
the present study of this quantity must be viewed as exploratory. However we note that in 
practice it is easy to tune L so that Eqn. El holds; one can thus expect to determine the N 
dependence of the mixing energy h in the near future. Theoretical arguments were given in 
section 2 why the mixing energy should be of order 1/A^ and therefore the dominant effect at 
large N in the cross-over region. 

The issue of mixing of near-degenerate states in the large expansion is an important 
and general one. For instance, a k = 2 string [2H1 of length L can be thought of as a state 
resulting from the mixing of two 'unperturbed' states^ which both have 1/A^^ corrections to 
their energy with respect to the planar limit value 2aL (a is the fundamental string tension). 
However a mixing energy of order 1/A^ implies that the lightest k = 2 string energy deviates 
by a term of that order from 2aL 38, . Thus studying the mixing strength of torelon pairs with 
glueballs in intermediate volumes satisfying Eg n . 1 1 71 addresses the same basic issue concerning 
SU(A^) gauge theories in a physical context where it is much easier for lattice calculations to 
give a clear-cut answer. 

To extend this work to L < 1/Tc, one has to confront the problem of long auto-correlation 
times on any quantity sensitive to the centre symmetry Z{N). This is especially true at 
large N, as outlined in the introduction. It then appears to be more advisable to use twisted 
boundary conditions, as was done for SU(2) in the work of Stephenson and Teper jSS] following 
the analytic work |37]. Working with larger gauge groups would then allow one to test the 
theoretical ideas of reduction. Eventually volumes are reached where perturbation theory 
becomes applicable and accurate (L ~ 0.2fm), so that the connexion with ab initio analytic 
methods jT| can be established. It would thus be interesting to investigate the spectrum 
of the effective small-volume Hamiltonian for the SU(A^) theory, A^ large. A necessary and 
non-trivial condition to give agreement with the large-volume spectrum, that is satisfied by 

^They correspond to operators Tr{P^} and {Tr{P})^, if P is the Polyakov loop. 
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this Hamiltonian even at finite A'^, is that it becomes spherically symmetric on very small 
volumes fTI. In view of the lessons learnt in this work and the fact that the multi-torelon 
operators decouple from the glueball-creating operators at large A^, we would expect only 
certain energy levels to match those of the gauge theory in infinite volume. 
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IR 


L = 8 


9 


10 


11 


12 


torelon 


0.2107(14) 


0.2599(30) 


0.3201(51) 


0.4016(48) 


0.4585(25) 


At+ 


0.54(1) 


0.53(1) 


0.55(1) 


0.63(1) 


0.661(9) 




0.731(8) 


0.77(2) 


0.77(2) 


0.88(3) 


0.93(2) 




1.23(4) 


1.24(7) 


1.15(3) 


< 1.37(2) 


1.24(5) 


E++ 


0.519(4) 


0.583(16) 


0.67(1) 


0.883(16) 


0.91(1) 


E++* 


1.12(2) 


1.07(3) 


1.04(1) 


0.97(3) 


1.075(7) 




1.39(4) 


1.46(9) 


1.42(4) 


1.54(12) 


1.46(7) 


rp+ + 

^2 


1.16(2) 


1.18(2) 


1.12(1) 


1.127(5) 


1.06(1) 


rp++* 


1.46(6) 


1.59(8) 


1.55(5) 


< 1.72(2) 


1.56(2) 




1.07(3) 


1.04(3) 


1.19(3) 


1.20(3) 


1.25(5) 


E~+ 


< 1.62(2) 


1.51(6) 


1.45(4) 


1.42(5) 


1.39(6) 




1.52(7) 


1.53(5) 


1.35(3) 


1.44(4) 


1.41(1) 


r+- 


1.50(6) 


1.46(5) 


1.41(3) 


< 1.48(3) 


1.39(3) 



IR 


L =14 


16 


20 


24 


torelon 0.5739(34) 


0.6781(33) 


0.8860(64) 


1.082(60) 




0.68(1) 


0.697(4) 


0.701(3) 


0.698(2) 


A++* 




1.17(2) 


< 1.32(3) 


< 1.23(2) 


E++ 


1.06(2) 


1.060(6) 


1.07(1) 


1.073(3) 


E++* 




1.43(1) 


< 1.56(4) 


< 1.48(2) 


rr+ + 

^2 




1.070(5) 






rn++* 




< 1.50(1) 










1.15(1) 






E-+ 




1.36(1) 










1.37(1) 








< 1.42(2) 


1.31(2) 


1.31(2) 


1.34(2) 



Table 1: The SU(3) spectrum (in lattice units) in finite volume from Monte-Carlo simulations 
at ^ = 6.0 (a ~ 0.097fm, V = x Lt). U is 24 for L > 12, 32 for L = 9 and 11, and 48 for 
L = 8 and 10. The method is described at the beginning of section |21 
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IR 


L = 8 


10 


16 


torelon 


0.2555(45) 


0.4514(40) 


0.8610(53) 




4859('69) 


6670f87) 


729fl0l 


^1 


0.7297(97) 


0.981(12) 


1.338(13) 


E++ 


5465f71) 


893fl4) 


1 139(11") 




1.1342(78) 


1.163(11) 


1.606(16) 




1.1682(68) 


1.152(10) 


1.143(16) 


rp++* 


1.370(30) 


1.599(22) 


1.622(16) 




0.924(14) 


1.060(72) 


0.991(86) 


E-+ 


1.528(18) 


1.503(18) 


1.507(50) 




1.525(37) 


1.477(40) 


1.423(41) 




1.491(20) 


1.439(21) 


1.453(15) 



Table 2: The spatial volume dependence of SU(4) glueball masses in lattice units at /3 = 10.90 
(a ~ O.llfm, 1/ = X U). U = 24 for L = 16, Lt = 36 for L = 8 and 10. 



IR 


^ 

L = 8 


10 


16 


20 


torelon 


0.3096(38) 


0.5016(63) 


0.930(10) 


1.2101(85) 




0.494(10) 


0.707(11) 


0.742(17) 


0.737(18) 




0.782(11) 


1.063(24) 


1.343(24) 


1.315(33) 


E++ 


0.589(12) 


0.942(27) 


1.153(18) 


1.180(12) 


E++* 


1.160(14) 


1.175(17) 






rn++ 


1.1786(88) 


1.177(12) 


1.167(13) 


1.160(10) 




0.941(55) 


1.203(26) 


<1. 127(51) 


<1.293(19) 


E-+ 


1.532(68) 


<1.573(28) 


1.498(21) 


<1.61(2) 




1.532(18) 


1.496(26) 


1.498(26) 


1.521(16) 




1.487(25) 


1.462(27) 


1.474(17) 


1.47(2) 



Table 3: The spatial volume dependence of SU(6) glueball masses in lattice units at /? = 25.05 
(a ~ O.llfm, V = L^ y. Lt). U = 24 for L = 16 and 20, U = 36 for L = 8 and 10. 
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IR 



L = 8 



10 



11 



12 



SU(3) 


. 


m+ 
rri- 


0.6394(44) 
0.1708(80) 


0.648(18) 
0.248(18) 


0.66(2) 
0.22(2) 


0.760(17) 
0.255(30) 


0.792(13) 
0.279(21) 




T— 1 1 1 


m+ 
rri- 


0.793(21) 
0.6182(93) 


0.794(37) 
0.478(28) 


0.855(14) 
0.37(1) 


0.918(35) 
0.120(20) 


0.9935(53) 
0.177(10) 


SU(4) 


. J L 


m+ 
rri-. 


0.601(11) 
0.237(8) 




0.8239(80) 
0.314(13) 








T— 1 J L 


m+ 
rri- 


0.8415(76) 
0.575(15) 




1.0268(75) 
0.263(12) 






SU(6) 


A++ 


m_ 


0.6379(81) 
0.285(12) 




0.885(14) 
0.355(24) 








E++ 


m_ 


0.8749(79) 
0.582(28) 




1.062(20) 
0.241(36) 







Table 4: The averages masses nij^ = {mQ+mi) /2 of and mass gaps m_ = mi — nio between the 
fundamental and first-excited states in the Af^ and E'^'^ sectors (in lattice units). Statistical 
correlations are taken into account in the standard way by the jackknife method for m+; for 
m_, we give the mass plateau of Ci{t)/Co{t), where Co{t), Ci{t) are the correlators of respec- 
tively the fundamental and first-excited operators obtained from the variational method |39j . 
Its variance is again determined by the jackknife method. 



Average masses in the A^"*""*" and E"*""*" sectors 




Figure 1: Average of the fundamental and first excited states in the Af^ and E^~^ sectors 
(black symbols). This energy level is independent of two-state mixing. Also plotted (white 
symbols) is mT{L) + mQ[L = oo). Where black and white symbol overlap, 'intrinsic' finite- 
volume effects are small. 
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Polyakov loop Monte-Carlo sequence V=16 x24 
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Polyakov loop Monte-Carlo sequence V=8 x48 
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Figure 2: Montc-Carlo sequence of the average value of the spatial-torelon operators over 10^ 
measurements; the measurements are separated by 4 sweeps. 
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Lightest 4D SU(3) Aj"*""*" states in finite volume 



h-1 
B 

03 




Lightest 4D SU(3) E"*""*" and T2^^ states in finite volume 



I 




Figure 3: Volume dependence of the lightest states in the representations that mix with torelon 
pairs. In the top plot, the smooth curve is a fit(0.699 — 674exp(— 0.605L)/L) discussed in the 
text. 
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4D SU(3) Aj"*" and Tj"*"' states in finite volume 
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4D SU(3) E"*" and T2"^ states in finite volume 




Figure 4: Volume dependence of some light states in various representations. 
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Figure 5: Finite- volume mass shifts m{L)/m(L = oo) — 1 in various representations; 
SU(3)— > black circles, SU(4)— white circles, SU(6)— white squares. The thin lines indicate 
the ±5% levels of variation. 
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